Decay rate and renormalized frequency shift of superradiant exciton in a cylindrical 

quantum wire 
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The decay rate and renormalized frequency shift of superradiant exciton in a cyhndrical quantum 
wire are studied. The transition behavior from ID wire to 2D film is examined through the property 
of the radiative decay. Similar to the case in a quantum well, the decay rate of the higher mode 
exciton is larger than that of the lower mode one. Moreover, it is also found the decay rate and 
frequency shift do not show oscillatory dependence on wire radius because of the conservation of 
angular momentum. 
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Since Dicke[|i| pointed out the concept of superradi- 
ance, the coherent effect for spontaneous radiation of 
various systems has attracted extensive interest both 
theoretically and experimentally ||, ||, |[ |j. In bulk 
crystal, the excitons will couple with photons to form 
polaritonsj^-mixed modes in which energy oscillates 
back and forth between the exciton and the radiation 
field. What makes the excitons trapped in the bulk crys- 
tal is the conservation of crystal momentum. If one con- 
siders a thin filniQ, the excitons can undergo radiative 
decay as a result of the broken crystal symmetry along 
the normal direction of the film plane. The decay rate of 
excitons in a thin film is enhanced by a factor of (A/d)^ 
compared to a lone exciton in an empty lattice, where A 
is the wave length of emitted photon and d is the lattice 
constant of the film. 

Lots of investigations on the radiative linewidth of ex- 
citons in quantum wells have been performed. An ab- 
normal increase of excitonic radiative lifetime with the 
decrease of well width below bnm for luajGai^xAs/InP 
quantum well was observed by Cebulla et al. . Brandt 
et al.^ measured the radiative lifetime of excitons in 
InAs quantum sheets and observed the increasing of ra- 
diative lifetime with the decreasing of well thickness. 
Hanamurajl^ investigated theoretically the radiative de- 
cay rate of quantum dot and quantum well. The obtained 
results are in agreement with that of Lee and Liu's 
prediction for thin films. Knoester|l^ studied the radia- 
tive dynamics crossover from the small thickness, super- 
radiant exciton regime to bulk crystal, polariton regime. 
The oscillating dependence of the radiative width of the 
exciton-like polaritons with the lowest energy on the crys- 
tal thickness was found. Recently, G. Bjork et aLfl^ 
examined the relationship between atomic and excitonic 
superradiance in thin and thick slab geometries. They 
demonstrated that superradiance can be treated by a uni- 
fied formalism for atoms, Frenkel excitons, and Wannier 
excitons. In V. M. Agranovich et al.'s work||l^, a de- 
tailed microscopic study of Frenkel exciton-polariton in 



crystal slabs of arbitrary thickness was performed. 

For lower dimensional systems, the decay rate of the 
exciton is enhanced by a factor of X/d in a linear 
chain iQ. First observation of superradiant short life- 
times of excitons was performed by Ya. Aaviksoo et 
g/. fisll on surface states of the anthracene crystal. A. L. 
Ivanov and H. Haug|l^ predicted the existence of exci- 
ton crystal, which favors coherent emission in the form of 
superradiance, in quantum wires. Y. Manabe et al.^ 
considered the superradiance of interacting Frenkel ex- 
citons in a linear chain. Recently, we have also shown 
the superradiant decay of the quantum wire exciton is 
greatly enhanced in a planar microcavityjl^. In this pa- 
per, the exciton is assumed to be confined in a hollow 
cylindrical quantum wire. The crossover of the superra- 
diant decay rate of the exciton from a wire with small 
radius to the 2D limit is explicitly obtained. Moreover, 
the crossover of the coherent frequency shift from narrow 
wire to thin film is also investigated by using the method 
of renormalization |l9| [20[ |. 

We consider a free exciton in a hollow cylindrical quan- 
tum wire lying on the z axis with simple cubic structure 
(Fig. 1). For the physical phenomenon we are interested 
in, we shall concentrate on the investigations of semicon- 
ductor quantum wires rather than existing carbon wires 
since the wavelength of the emitted photon is 



ni 22 



'Electronic address: 
t Electronic address: 



ynchen.eD87g@nctu.cdu.tw 



ischuu@cc.nctu.cdu.tw 



much larger than the diameter of the single-wall carbon 
tubes. The crossover behavior from ID wire to 2D film 
may not be examined on carbon systems. However, we 
hope our model can serve as a first step toward the un- 
derstanding of the exciton decay in carbon nanotubes. 

If the difference between the inner radius (p^) and 
outer radius (p>) is much smaller than the effective Bohr 
radius of the exciton, i. e. AL << Uex, one may ap- 
proximate that the exciton is confined on the cylindrical 
surface with radius p(w w P>)- This means the ex- 
citon is trapped in an infinite deep and narrow quantum 
wire. The radius p is about Nd/2Tr, where N is the num- 
ber of lattice points in the circumference direction and d 
is the lattice spacing. If one further assumes the radius of 
the wire is much larger than the effective Bohr radius of 
the exciton, variations of the wire radius only cause few 
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where m is the electron mass, {p^,ip^,li) is the position 
of the electron i in the cylindrical wire, Pi is the corre- 
sponding momentum operator of the electron i , eq'^'k'_ 

is the polarization vector of the photon, and i?^)^' (q'p) is 
the Hankel function. 

The essential quantity involved is the matrix element 
of H' between the ground state |G) and the exciton state 
\i^,kz). We know that the interaction matrix elements of 
H' can be written as 



FIG. 1: Schematic view of the quantum wire structure and 
its defining potential profile. 
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because the exciton state can be expressed as 



(4) 



changes on the Wannier exciton wavefunction. In this 
case, the main contributions to the superradiant decay 
rate and frequency shift still come from the bandgap en- 
ergy and the number of the lattice points within a wave- 
length of the emitted photon[^. Therefore, one can first 
consider the exciton as a particle with angular momen- 
tum V and longitudinal momentum kz- After figuring 
out the decay rate and frequency shift of a Frenkel exci- 
ton, the corresponding ones of a Wannier Exciton can be 
obtained by replacing the single-atom dipole matrix ele- 
ment X with the effective dipole matrix element [Q. Thus, 
the Hamiltonian for the exciton is 



in which the excited state |c, (Z, ip); v, {I, (p)) is defined as 



|c, {I, p);v, {I, Lp)) = a);^(,^^)a„,(;,<p) \G) 



(6) 



where a)^ ^■^iO'v,(i,>fi)) the creation (destruction) oper- 
ator of an electron in the conduction (valence) band at 
lattice site {1,^). The expansion coefficient J7*j, (Ij'p) is 
the exciton wave function in the cylindrical tubule: 
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where cj,j, and c^k^ are the creation and destruction op- 
erators of the exciton, respectively. The Hamiltonian of 
free photon is 



k?Y'%>.>k'S''^'k'^, (2) 



q' l/'k'^ 



where b^q,^,/,, and 6gv'fc^ are, respectively, the creation 
and destruction operators of the photon. The wave vector 
k' of the photon is separated into two parts: k'^ is the 
parallel component of k' on the z direction such that 

The interaction between the exciton and the photon 
can be expressed as 



u:kAi,v) = 



1 I 
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where the coefficient l/\/]V' is for the normalization of 
the state fc^). 

After summing over I and ip in Eq. (4), we have 
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where 



A^+„^,fc,+ff = V NN' J (Pt exp(z(/c2 + g)Tz + + n^) 

XWc{Tz,T^){-ih'V)w^{Tz,T^), 



FIG. 2: Decay Rate of the superradiant exciton as a function 
of the wire radius. The vertical and horizontal units are ps~^ 
and nm, respectively. And the solid, dashed, and dot-dashed 
lines represnet the f = 0,1, and 2 modes, respectively. In this 
and following graph the paramaters of a GaAs quantum well 
are chosen to obtain the numerical results. 

In Fig. 2 wc plot the decay rate j^^. as a function 
of wire radius. Our numerical results are obtained by 
employing the data of GaAs (band gap is 1.52 eV). For 
1/ = mode(solid line), the decay rate increases with the 
increasing of wire radius and approaches 2D limit as p is 
large comparing to the wavelength of the emitted photon 
in a 2D thin film. For higher modes(z^ = 1, 2), the results 
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depend on the kinetic energy '-^'^ ^^^"^ from 

the figure, the decay rate also approaches 2D hmit for 
large radius. On the other hand, the decay rate increases 
as the wire radius is decreased until the maximum decay 
rate is attained. After that a further reducing of the ra- 
dius leads to a sharp decrease of the decay rate. This is 
similar to the transition from 2D to 3D: the higher wave 
number modes have larger maximum decay rate[p^. One 
also notes in small radius regime, the decay rate of the 
higher mode is smaller than that of the lower mode, i.e. 
lv=2,k. > lu=i,k. > lu=o,k,- This phenomenon is also 
found in the quantum well system, and can be ascribed to 
the effects of interference |2^. Another interesting result 
is that the decay rate does not show any oscillatory de- 
pendence on radius while it shows oscillatory dependence 
on layer thickness in the case of the semiconductor thin 
filmpl|, The reason is that the angular momentum 
of the cxciton is conserved in a cylindrical system, but, 
in a semiconductor thin film, the momentum is not con- 
served in the direction of broken symmetry. If we sum 
over all the sites in Eq. (4), the phase will run through 
a circle and thus preserve the same value. 

One might suspect the practical value of our numeri- 
cal results because our model is based on the one-layer 
assumption and the idealized assumptions , e.g. perfect 
cylindrical symmetry and infinite confinement. However, 
as we mentioned above, for the superradiant decay rate 
and frequency shift one only needs to replace the dipole 
matrix element with an effective one. This is why our 
result agrees well with that of a realistic GaAs quantum 
well in the large radius limit ||l0|| . On the other hand, 
such an assumption may slightly deviate from the real 
case as the radius is small (e.g. as the radius approaches 
ttex)- This is because the wavefunction of the exciton 
becomes more compact in the quantum wire with small 
radius limit, and it causes the increasing of the dipole 
matrix element x in Eq. (17). Besides, the quality of 
the quantum wire also influences the observation of the 
crossover, i.e. the thickness fluctuations should be con- 
trolled well, otherwise it will destroy the coherence in the 
circular direction. 

Now let us turn to the results for the renormalized 
frequency shift. The frequency shift in Eq. (15) can be 
expressed as 
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As seen from above, the frequency shift suffers from ul- 
traviolet divergence when g and are large, and has in- 
frared divergence when the denominator approaches zero. 
Following the procedure as shown in the work of Lee, 
Chuu, and Mei||l^, the divergent problem is solved by 
renormalization. 



In the usual renormalization procedure used in quan- 
tum field theory, the infinite quantities in Green function 
is subtracted by some infinite quantities to make them 
finite. The subtraction procedure is possible by sub- 
stituting some finite physical quantities, such as renor- 
malized masses, charges, and wave functions. The finite 
physical quantities observed from bare infinite quantities 
can be visualized, in the case of bare charge for exam- 
ple, as the polarization of vacuum. This is equivalent 
to say that there is virtual photon and opposite charged 
electron cloud surrounding the bare charge making the 
charge measured from outside become finite~a phenom- 
ena similar to that of shielding in dielectric material. But 
the renormalization procedure in quantum field theory is 
for free electrons. In the case of the condensed mat- 
ter, the electrons are confined by periodic potential and 
complex interactions. Borrowing from the concepts of 
renormalization used in quantum field theory, we have 



kf) — >0,a^oo 



(20) 



where the two limiting processes fco ^ and d ^ oo 
reduce the exciton to a free electron. In the limiting pro- 
cess d —^ oo, the ordinary exciton becomes a lone exciton 
standing alone in an empty lattice with no interaction 
with other atoms. And the process fco ^ means that 
there is no energy difference between electron and hole. 

We will now show that the ultraviolet divergence comes 
from the inclusion of Umklapp process. Define 



f^i^fc, (^¥')= J„kAk2,n^)exp{ik2l + in^(p) (21) 
with 



JvkAk2,n^ 
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From the above equations, we have 



^vk, {I, v) exp(-ikj - ivif) 
J2 J-^kAk2,n^) J2 e'('^^-'^^)'+''("--'')'^ 



(22) 



k^n^ 



gn^ 



k^ ■ 



(23) 



And from Eq. (21) we have lim^j^o^i/fc^ = ^i^k^ {I = 
0, (/J = 0), so we can rewrite the renormalized frequency 
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shift nil'l as 



with 



and 



hm f7^fe,(0,0) 

fco— >0 



(24) 



(25) 



(26) 



As can be seen from Eq. (25), the renormalization af- 
fects only the part Slj^fe^ (0, 0)-the frequency shift of the 
lone exciton in an empty lattice-and thus nothing to do 
with the correlation within the crystals. The coherent 
part VL""^] in Eqs. (24) and (26) is not touched by renor- 
malization procedure. The separation of ^uk^ into two 
parts as shown in Eq. (24) is conceptually equivalent 
to singling out of the source term of quantum electrody- 
namical divergence in a correlated system. 

Now we will investigate the origin of ultraviolet di- 
vergence. From Eq. (19) with the substitution — s- 

/(^andE,-/7^, wehave 



(27r/L, 



Airm^c^ 4-7 ^ ^ J " Jn'^ + (fc, + g)2 



V?'' + (fc. + 9? 



H, 



(1) 



(27) 



It is noted that the result in Eq. (27) is the same as 
^^i'fc^(0,0) from Eqs. (21) and (23). Hence we conclude 
that the ultraviolet divergence of fiyfe^ (0, 0) is really the 
same as that from Umklapp process of large g and Ui, that 
contribute to the full fii/fc^. Once £7^^,(0,0) of the lone 
exciton is rendered finite by renormalization via Eq. (25) , 
the Umklapp processes of large g and n that arise from 
the unrestricted sum over / and ip will also be rendered 
finite simultaneously. Accordingly, fi^^" can be written 



(28) 

and thus the ultraviolet divergence problem is solved. 

In the kz'-^O and v = mode, the renormalized result 
can be reduced as 
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As seen from above, the frequency shift suffers from 
infrared divergence when g' ^ or ~ /cq- This can be 
overcome by substituting —iTiW by —imcq'T(Rei. 19) in 
Eq. (18) when q' is small. It is equivalent to the dipole 
interaction form, if' ~ r • E. With this treatment, we 
have 



k^ ■ 



Hl^\q'p) 



ko - q' 



(30) 



with 
where 



7re h 



|e<?vfe, • x\ , when q' is large ^^^^ 
\^q'vk^ ■ , when q' is small 



(32) 



Eq. (30) can not be evaluated analytically. But for 
large radius, the asymptotic form of Hankel function is 
: H'^\q' p) ~ ■y/2/(7r(j''p)e*'' And the renormalized fre- 
quency shift reduces to 2D limit p9|: 



= -7single(r;^) 



(33) 



where 



7s 



2e2 E^, 



he 



k^d\ 



(34) 



is the radiative decay rate of a single isolated exciton. 
Analogous to the decay rate, the renormalized frequency 
shift is explicitly seen to be coherently enhanced by the 
same factor (l/fcod)^ as a result of the interaction of 
the phase-matched photon amplitude with the delocal- 
ized excitonic amplitude in the plane. In Fig. 3, we 
numerically calculated the frequency shift as a function 
of wire radius. One can see from the figure the renormal- 
ized frequency shift does not show oscillatory dependence 
on radius, either. For v = Q mode, unlike the behavior 
of decay rate, the magnitude of the frequency shift first 
increases with the decreasing of radius. After reaching a 
minimum point, the frequency shift approaches to zero 
rapidly. 

For usual semiconductors, the superradiant enhanced 
factor is about 10^ for excitons in the optical range. How- 
ever, due to the extreme smallness of 7singie itself, obser- 
vation of ^v,k^ is not expected to be easy. Its depen- 
dence on wire radius may be a useful feature to observe 
this quantity. Recently, R. A. Romer and M. E. Raikh 
studied theoretically the exciton absorption shredded by 



FIG. 3: Renormalized frequency shift of the superradiant ex- 
citon as a function of the wire radius. The vertical and hori- 
zontal units are ps^^ and nm, respectively. 

a magnetic flux $ in a quantum ring[^. They found 
the oscillator strength of the exciton is most enhanced 
when $ is equal to half of the universal flux quantum 
$0 = hc/e. And the oscillation period is equal to $0. 
Owing to the similarity in cylindrical geometry, the de- 



cay rate and the frequency shift may be observed ex- 
perimentally if one varies the magnetic flux or the wire 
radius. 

In summary, we have calculated the decay rate of the 
superradiant exciton in a hollow cylindrical quantum 
wire. Similar to the case in a quantum well, the higher 
wave number modes are shown to have larger maximum 
decay rate. It is also found the decay rate does not show 
any oscillatory dependence on wire radius because of the 
conservation of the angular momentum in the cylindrical 
system. On the other hand , the frequency shift of the 
exciton is properly renormalized by removing the ultra- 
violet and infrared divergence. It is found the shift does 
not shown oscillatory dependence on the radius, either. 
Some distinguishing features are pointed out and may be 
observed in a suitably designed experiment. 

We would like to thank Professor M. F. Lin of Na- 
tional Cheng Kung University for a helpful discussion. 
This work is supported partially by the National Science 
Council, Taiwan under the grant number NSC 91-2112- 
M-009-012. 



[1] R. H. Dicke, Phys. Rev. 93, 99 (1954). 

[2] N. Skiribanowitz, I. P. Herman, J. C. MacGillivrary, and 

M. S. Feld, Phys. Rev. Lett. 30, 309 (1973). 
[3] M. Gross, C. Fabre, P. Pillet, and S. Haroche, Phys. Rev. 

Lett. 36, 1035 (1976). 
[4] V. Ernst and P. Stehle, Phys. Rev. 176, 1456 (1968). 
[5] Y. C. Lee and D. L. Lin, Phys. Rev. 183, 147 (1969); 

183, 150 (1969). 
[6] J. J. Hopfield, Phys. Rev. 112, 1555 (1958). 
[7] K. C. Liu and Y. C. Lee, Physica 102A, 131 (1980). 
[8] U. CebuUa, G. Bacher, A. Forchel, G. Mayer, and W. T. 

Tsang, Phys. Rev. B 39, 6257 (1989). 
[9] O. Brandt, G. C. La Rocca, A. Heberle, A. Ruiz, and K. 
Ploog, Phys. Rev. B 45, 3803 (1992). 
[10] E. Hanamura, Phys. Rev. B 38, 1228 (1988). 
[11] J. Knoester, Phys. Rev. Lett. 68, 654 (1992). 
[12] G. Bjork , S. Pau, J. M. Jacobson, H. Cao, and Y. Ya- 

mamoto, Phys. Rev. B 52, 17310 (1995). 
[13] V. M. Agranovich, D. M. Basko, and O. A. Dubovsky, J. 

Chem. Phys. 106, 3896 (1997). 
[14] V. Agranovich and O. Dubovsky, JETP Letters 3, 223 



(1966). 

[15] Ya. Aaviksoo, Ya. Lippmaa, and T. Reinot, Optics and 

Spectroscopy (USSR) 62, 419 (1987). 
[16] A. L. Ivanov and H. Haug, Phys. Rev. Lett. 71, 3182 

(1993). 

[17] Y. Manabe, T. Tokihiro, and E. Hanamura, Phys. Rev. 

B 48, 2773 (1993). 
[18] Y. N. Chen and D. S. Chuu, T. Brandes, and B. Kramer, 

Phys. Rev. B 64, 125307 (2001). 
[19] Y. C. Lee, D. S. Chuu, and W. N. Mei, Phys. Rev. Lett 

69, 1081 (1992). 
[20] Y. N. Chen and D. S. Chuu, Phys. Rev. B 61, 10815 

(2000). 

[21] K. Suenaga, C. Colliex , N. Demoncy, A. Loiseau,H. Pas- 
card, and F. Willaime, Science 278, 653 (1997). 

[22] Y. Zhang, K. Suenaga, C. Colliex, and S. lijima. Science 
281, 973 (1998). 

[23] R. A. Romer and M. E. Raikh, Phys. Rev. B 62, 7045 
(2000). 



